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The theory of the preroughening transition of an unreconstructed surface, and the ensuing 
disordered fiat (DOF) phase, is formulated in terms of interacting steps. Finite terraces play a 
crucial role in the formulation. We start by mapping the statistical mechanics of interacting (up 
and down) steps onto the quantum mechanics of two species of one-dimensional hard-core bosons. 
The effect of finite terraces translates into a number-non-conserving term in the boson Hamiltonian, 
which does not allow a description in terms of fermions, but leads to a two-chain spin problem. The 
Heisenberg spin-1 chain is recovered as a special limiting case. The global phase diagram is rich. 
We find the DOF phase is stabilized by short-range repulsions of like steps. On-site repulsion of 
up-down steps is essential in producing a DOF phase, whereas an off-site attraction between them 
is favorable but not required. Step-step correlation functions and terrace width distributions can be 
directly calculated with this method. 
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I. INTRODUCTION 



The surface roughening transition and the nature of the rough phase are theoretically very well understood by a 
variety of approaches ranging from phcnomcnological descriptions, based on the sine-Gordon model, to microscopic 
solid-on-solid (SOS) models, g] 

The preroughening transition (PR) and the ensuing disordered flat phase (DOF), both predicted several years 
ago by Rommelse and den Nijs, |2| have also been studied and characterized within certain restricted solid-on-solid 
(RSOS) models. ||-|!| Although the physics behind these models, and the ingredients stabilizing the DOF phase, have 
been discussed in some detail, it is still useful to explore this subject from a different, and perhaps more physically 
appealing, perspective. The RSOS models, in particular, do not directly emphasize steps, terraces, and kinks, which on 
the other hand are very crucial actors in these transitions. This need is made more urgent by the recent experimental 
evidences for preroughening on rare gas solid (111) surfaces, J^Jll[] which calls for a detailed reinvestigation of the 
step-step interactions, M or the reconstructive tendencies, crucial to obtain a DOF phase. 



In the context of roughening, Villain and Vilfan, 12]-14|] den Nijs, p"5| and Balents and Kardar, |l(| for the case of 
reconstructed surfaces, and Villain, Grempel, and Lapujoulade, |17| for the case of vicinal surfaces, have shown that 
a more phcnomcnological approach based on working directly with steps yields a very direct picture of the physics 
involved. 

Anisotropic surfaces, in particular, have a definite direction of stronger bonding along which a step tends to run; 
kinks on such steps, involving the breaking of strong bonds, are energetically expensive. The (110) face of fee noble 
metals (Au,Pt,Ag,. . .) is a physical realization of an anisotropic surface, with a stronger, compact, [110] direction, and 
a softer [001] direction. In the strongly anisotropic limit, the transfer matrix problem for the system of steps can be 
mapped onto the imaginary-time evolution of a system of quantum particles in one dimension (ID), the imaginary- 
time being the preferred direction in which the steps run. ]l8|-p0| This is a well known mapping, hea vily exploited, 
for instance, in the theory of uniaxial commensurate-to-incommensurate transitions of adsorbates. |20| , |l9| 

Following a line of thought initiated by Ref. |l3), Balents and Kardar applied the full machinery of the theory 
of interacting fermions in one dimension to explore the possible phase diagrams of generic (px 1) reconstructed 
anisotropic surfaces. |ll| In their approach, steps of double height are forbidden (as energetically too expensive), 
while up and down monoatomic steps are mapped onto spin-1/2 fermions in ID, described (in the continuum limit) 
by the hamiltonian fllq] 



H = 



E 



+ / dx dx' 'n a (x)V a a' (x — x')n a i (x') 



(1) 



Here 7 is the inverse line tension of a step, V aa > (x) is an interaction between steps, and the remaining notation 
is standard. This Hamiltonian describes infinite steps, traversing the entire length of the sample, as the number 
of particles is conserved by the (imaginary-time) dynamics. In reality, steps on surfaces can lead to finite defects 
by forming loops (i.e., finite terraces) on the surface. ||l^ , |l6f The order of the reconstruction p, dictates, through 
the symmetry of the different ground states, the form of the "loop" terms which are allowed. jl3];[l6) For a {p x 1) 
reconstructed surface, the Hamiltonian H has to be supplemented with a term of the type 



-Sloop = X dx 



■p-i 



Y\_tp^(x + ka)tpi(x + ka) + H.c. 



,fc=0 



(2) 



where a is a cut-off distance of the order of the lattice constant. Balents and Kardar argued, by power counting, 
that the effect of finite terraces, i.e., the introduction of -Sloop, is irrelevant for p > 2, marginal for p = 2, and 
strongly relevant for p = 1. They went on by addressing in detail the p = 2 case, of relevance to the (110) missing-row 
reconstructed facet of Au. (See also Ref. for closely related work on the p — 2 case.) However, the unreconstructed 
(p = 1) case was not pursued further. 

The approach we take in the present paper is similar in spirit. Our specific goal, however, is to address the 
question of the presence of a DOF phase, and the preroughening transition, for unreconstructed surfaces. Thus, in 
the classification introduced above, we are now interested in detail in the p — 1 case. Technically this leads, as we 
shall see, to significant differences with respect, for instance, to Ref. p6[ , and to a new phase diagram quite different 
from the p > 1 cases. 

The crucial point is that for p = 1 the loop terms (see Eq. (Q)) are of the BCS-like form A J dx(ip^ (%)ipl i x ) + H.c), 
i.e., a strongly relevant one-body piece. This might appear as just a minor complication, at first glance, since quadratic 
terms can be easily diagonalized by a Bogoliubov transformation. Closer consideration, however, leads to reconsider 
the whole mapping. Fermionic minus signs have no role whatsoever in the mapping of a classical statistical mechanics 
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problem. The natural statistics to use is always the bosonic. In the present case, a hard-core constraint will be 
necessary, in order to implement the appropriate configurational space (for instance, the non-crossing constraint 
for steps of the same type). In a one-dimensional quantum problem, the choice of the statistics is, quite often, 
not a big problem, as we can transform, by a Wigner-Jordan transformation, hard-core bosons into fermions, with a 
transformed hamiltonian which has exactly the same form (only boundary conditions have to be considered carefully) . 
In our case, however, pairing terms of the type a^-fffli+i 4, which (see below) are essential to describe finite terraces, 
do not transform into simple fermionic BCS-like terms, and become non-local after a Wigner-Jordan transformation. 
This will force us to work with hard-core bosons. 

Our approach, in summary, is as follows. We assume, as in Ref. ||l6), that the only relevant extended defects are 
monoatomic steps, which can be either up or down. Steps of the same kind are forbidden to cross, while steps of 
different type can cross. Moreover, steps interact with each other, have kinks, and can form finite terraces on the 
surface. These steps are then mapped onto world-lines of hard-core bosons in one dimension. Kinks on the steps 
correspond to hopping terms in the quantum Hamiltonian. Pairs of up-down steps which are created and annihilated 
to form finite terraces on the surface, give number-non-conserving terms in the quantum Hamiltonian. (2^] Pairwise 
interactions between the steps are taken into account by corresponding two-body terms in the quantum model. 

The present work is concerned with the case of a low-index unreconstructed surface. Extensions to the case of 
vicinals, for which the long-ranged nature of the step-step interactions is an essential ingredient, are left to a future 
study. 

Our main goals in working out this type of approach to PR are the following: (a) to build a formulation providing 
a more direct access to the physics of PR, which is somewhat hidden in the RSOS formulations; (b) to explore more 
directly the role of step-step interactions; (c) to study step-step correlation functions and terrace width distributions, 
not available so far. As it turns out, we have found that this approach is quite successful on all three accounts. 

The paper is organized as follows. In Sec. || we present in detail the classical statistical mechanics model of 



interacting steps, which is then mapped onto the corresponding ID quantum model of hard-core bosons in Sec. [II. In 



Sec. IV we consider in detail the spin-1 limit, obtained by setting the on-site step-step repulsion to infinity, and then 
map the general case to a problem of two spin-1/2 Heisenberg chains. Sec. contains a summary of bosonization 



plus finite-size scaling calculations done in order to extract the phase diag ram of the model. Sec. VI summarizes the 



relevant order parameters and correlation functions investigated. In Sec. |VII| we present our results for the overall 



phase diagram of the model. Sec. VIIl illustrates our results for the step-step correlations and the terrace width 



distributions. Finally, Sec. IX contains a discussion of the results and some conclusive remarks. 



II. THE MODEL 

We assume the only relevant extended defects involved in the surface PR transition to be steps, which can be either 
up or down. These steps interact with each other, they can have kinks, and they can form finite terraces on the 
surface. Fig. |I| shows a schematic picture of a surface with steps. 

Our model will be defined on a square lattice, and we will assume the steps to run preferentially in one direction (the 
vertical direction in Fig. |l|). The steps are only allowed to make simple nearest neighbor kinks. Steps running in the 
horizontal direction are assumed to be energetically expensive and neglected. Hence, our surface is, by construction, 
very highly anisotropic. 

We define the model by assigning its transfer matrix along y. Denoting by |S(j)) and |S(j + 1)) the configuration of 
the j-th and j + 1-th horizontal strips, we have 

(S(j + 1)|T |S(j)> = 

= exp (-0 {isN® + 8 K N ( ^ +1) + £ + 5 ex N^ 0+1) + W*»p) j (3) 

where (3 is the inverse temperature, and (see also Fig. ||): 

• 5s represents the energy cost (per unit length) of a step running along the y-direction; N s ^ is the number of 
steps in the strip j; 

• SkN^'^ 1 ^ is the energy cost of N^ +1 ' kinks between strip j and j + 1; 

• 5t s -^t ' J+1 *' is the energy cost for the creation of N^ +1 ' terraces of "size" s between strip j and j + 1. We will 
always assume s = 1, or (see Fig. 0). 
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• 5 ex N^ +1 ^ is the energy associated to the crossing of N^ +1 ^ pairs of opposite steps between strip j and j + 1; 

• Kitep-stcp — ^" + V ± > with V"ll and V 1 - describing respectively the interaction between steps of the same kind 
and of the opposite kind. For V", we assume a generic repulsive interaction 

K " = Y Y Vk-i n ^n k , a , (4) 

a k>i 

with VjjLj possibly possessing an elastic long-range tail of the form \k — i\~ 2 . Here Jfy.-m) is 1 if there is a step 
up (down) at site i. Similarly, we assume V 1 " to be given by 

y± = v o~ Y n t,l n t,l + Y Y Vk-in iltT n k> a • (5) 

i (7 k>i 

The sign of the terms in V , particularly at short-range, depends on microscopic details and need not be 
specified at this stage. 

If we assume periodic boundary conditions in y direction, i.e., |S(N y + 1)) = |S(1)), the partition function of this 
system is 

Z = lim Tr T N * . (6) 

N y — »00 



III. THE QUANTUM MODEL 



It is well known that, in the strong anisotropy (or time-continuum) limit, many D dimensional classical problems 
can be mapped onto D — 1 dimensional quantum problems. |Tg| ] This relationship is established by means of the path 
integral formalism. In particular, the up and down steps are mathematically equivalent to the world lines of spin-up 
and spin-down hard-core bosons, and the preferential direction in which the steps run plays the role of time in the 
quantum problem. The hard-core condition is imposed in order to implement the non-crossing condition for steps 
of the same type, a physically justified restriction, in view of the large energetic cost of double-step regions. The 
non-crossing constraint for steps of the same type would be automatically satisfied by the Pauli principle if we were 
to deal with spin-1/2 fermions. (See below for more comments on the problem of quantum statistics.) 

We consider the following quantum Hamiltonian 

H = -t 11 ^(4 CT a i+ i, CT + H.c.) - fiN + 



1 a f 



+ H.c.)-f 



i+1,5 



H.c.) -t e 



Y^ 



+ YY V }-i n i>° n i,° + V Y n *T"U + Y Y V 5~-i n i><r n 3,° ( ? ) 

with di^a representing the destruction operator for a spin a hard-core boson and with N = + iVj_ the total number 
of particles. We will work in the subspace N-\ = Ni, for a low- index surface. (For a vicinal surface of angle </>, we 
would have (N± — N-\) — Ltan</>.) 

Within a path-integral approach, it can be shown Jl8| that the ground state properties of this quantum Hamiltonian 
correspond to the temperature properties of the classical step model, whose transfer matrix is given by Eq. (^), in the 
large anisotropy limit. Specifically, the classical parameters turn out to be given by: 



-0S K 



et 



0,1 ' 
(11,-L) 



et e 



11,-L) 



(8) 
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where e is the Trotter discretization time for the quantum path-integral. 

The mapping is asymptotically correct only in the limit e — > 0. This a) implies, clearly, a strong anisotropy limit for 
the classical problem and b) does not allow a straightforward identification of a classical low- or high-temperature 
limit. Indeed, if all the parameters of the quantum problem are of order one, taking f3 — ► or (3 —> oo makes Eq. (^) 
incompatible with the requirement that the left-hand side should be a small quantity (of order e). In other words, 
the mapping is justified so long as 

Kinetic couplings = (Sk, <5t , • • •) ^> T ^> (5s, V, ■ ■ •) = Potential couplings , 

and nothing can be said, in principle, about the infinite temperature limit. This should be always kept in mind when 



considering the infinite temperature limit from the quantum model point of view (see, e.g., sec. VII B). 

It is also worth stressing that the (hard-core) boson statistics of the a operators in the Hamiltonian is crucial to 
the nature of the phases and transition lines in the phase diagram. Indeed, unlike other terms in the Hamiltonian, 
the terrace creation terms cannot be translated into simple (i.e., local) fermionic BCS-like terms by a Jordan- Wigner 
transformation. In such instances, the correct statistics to use is undoubtedly the bosonic one, as fermionic minus 
signs do not appear in a classical statistical mechanics problem. This point seems to be not always appreciated in the 
literature. §f] 

We now make contact with previous work in the context of surface physics. The model in Eq. (pi), with t ex = and 
^i = 0' nas been considered, in its continuum version, by Balents and Kardar. (See also Rcfs. Q and jl5| for 
related work.) In the absence of t 1 - -terms, particles are taken to be fermions. The emphasis of Ref. |l6| was on (p x 1) 
reconstructed surfaces, particularly with p = 2. The effect of finite terraces, i. enclosed loops of steps, was argued to 
be irrelevant for p > 2, marginal for p = 2, and strongly relevant for p — 1. [|16| The unreconstructed (p = 1) case, 
however, was not pursued at all. As just argued, the p = 1 case cannot be tackled in terms of fermions. The effect of 
the finite terraces on an unreconstructed surface - the t 1 - terms in Eq. (Q) - is one of the points addressed in detail in 
the present work. Moreover, we show that restricting the analysis to a simple Hubbard-type on-site interaction does 
not lead to the full richness of the phase diagram; nearest-neighbor interactions are essential in order to stabilize, for 
instance, a DOF phase. 



IV. CANONICAL TRANSFORMATIONS 

A. Vq 1 = oo: Mapping to a spin-1 chain. 

For a special choice of the parameters, the model in Eq. (Q) reduces to a well studied problem. Consider the case in 
which both V" and are truncated to nearest neighbors, vj\^ — v\ 5j t i+x, in the limit of infinite on-site repulsion 
of opposite steps, — > oo. The limit Vq — > oo enforces, in absence of t ex , a non-crossing condition for opposite 
steps as well, and allows only three states per site, which we can easily map onto a spin-1 variable as follows: 

step up at T |0> = \S? = +1) 
no step |0) = |5f = 0) 
step down at JO) = |5f = -1) . (9) 

It is then straightforward to show that all possible matrix elements of H in Eq. (0) coincide exactly with those of the 
spin-1 Heisenberg chain 

#Hcis = -^f YK S tSi +1 + H.c.) + J Z J2 S*St +1 + D Y,(SD 2 , (10) 

i i i 

if the following parameter choice is made: 

i 11 =ti = J xy , t ex = , V} = -Vf = J Z , -11 = D . (11) 
A nonzero t ex term, when present, translates into a quartic spin term 

-t ex Y l [(S+nS7 +1 f + K.c.}, (12) 
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and will be shown to be relevant in stabilizing the rough phase for finite repulsive V± . The phase diagram for this 
special case, in the surface physics relevant region — fx = D > 0, can be directly borrowed from the literature. p|,pH 

The Heisenberg spin-1 chain was also obtained by another route, namely by the quantum mapping of RSOS models, 
by den Nijs and Rommelse, who gave a very detailed discussion of the surface physics interpretation of the different 
phases. || The region in which a DOF phase is stabilized is found for J z > 0, i.e., it corresponds to repulsive and 
attractive (see Eq. (pi])). The latter condition is, however, not crucial, as we shall see later on. The flatness of 
the DOF phase is directly related to the Haldane gap in the spin-1 chain. |3| If the interactions are much bigger than 
the cost of a unit of step (J z S> D, see line (a) in Fig. |3|), at low temperatures (i.e., for large J z /J xy ) the system will 
reconstruct into an ordered sequence of up-down steps (the Necl phase of the spin-1 chain). By increasing temperature 
(i.e., lowering J z /J xy ) the system undergoes an Ising transition to a DOF phase, in which the positional order of the 
up-down sequence of steps is lost. || If the cost of the unit step is larger than the interactions, ( J z <C D, see line (b) 
in Fig. ||), the low temperature phase is flat, and the (preroughening) transition to the DOF phase has non- universal 
exponents. J^] 

If we impose the condition — ► oo without assuming Eq. (|ll|), what we are considering is always a three-state- 
per-site problem, but the resulting Hamiltonian does not have the simple bilinear form ( |l0|) in terms of the spin-1 
operators. We will refer to this general case as a spin-1 chain. The specific case in Eq. ( pX)[ ) will be referred to as 
Heisenberg spin-1 chain. 



B. The general case: Mapping to two spin-1/2 chains 

In order to study the model in Eq. (|t]) for more general parameter values, in particular for ^ oo, it is convenient 
to abandon the spin-1 representation and map onto a problem of two coupled spin-1/2 chains, where the total number 
of states per rung is four, instead of three. Introducing the usual spin-1/2 representation for each of the two species 
of hard-core bosons, 



a 



ltd) ~~ ^1(2) ' n i,t(i) - ^,1(2) + 2 



and performing a 7r-rotation for the spins around the x-axis, which amounts to a particle-hole transformation for 
the down-bosons, S^ 2 ~ * Sf 2 an d Sf 2 — + ~Sf 2 , one can rewrite the Hamiltonian (0) as the following model of two 
spin-1/2 chains (a = 1,2 denoting the chain) with opposite magnetic fields: 

H = -*" E(Oi~+l,a + H.C.) + V} J2 Sl a SUl,a + hY,( S tl " S h) 

i,a i,a i 

-to E( 5 i^2 + H - c -) - V E s ii s h + 

i i 

-*t E ls v,a.i,-. + H - c -) - v ^ E s ust + iA + 

— lex E (^+1,1^1,1 ^1+ 1,2 ^4,2 + H.C.) , (13) 

i 

where, for simplicity, we have considered only interaction terms up to nearest neighbors and we have omitted the 
constant ^[V rX (g = 0) + V"(q = 0)] (with V^ ±s> (q) we denote the Fourier transforms of the potentials, and V"^(g = 

0) = J2i Vj-i )■ The magnetic field h is related to the chemical potential fi in the following way: 



V X (q = 0) + VH (q = 0) =-n+ + V x + V?) 



Notice that after the spin-rotation, performed to get a standard S 2 coupling between the chains starting from the 
ajaj boson term, the signs of the V L terms are all changed, since Sf >2 — > —Sf 2 . For the same reason, the chemical 
potential terms transform into opposite magnetic fields for the two chains. The condition iVj = N± reads, after the 
canonical transformation, as zero total magnetization for the spins, ^ (S^i + $i2) = 0- P3i 



G 



V. LOW-ENERGY HAMILTONIAN FROM FINITE-SIZE DIAGONALIZATION 



With a well defined quantum spin chain problem at hand, it is standard practice to study its phase diagram by a 
combination of field-theoretical arguments and finite-size exact diagonalization data. 

At weak coupling, a standard field-theory approach to one dimensional quantum systems consists in applying 
bosonization techniques. This was done by Strong and Millis for the case of two coupled spin-1 /2 Heisenberg chains, 
i.e., a special case of Eq. (13) with t± = = 0, t ex = 0, h = 0, and Vq~ = 2<q = — Jk- [p3| The procedure can 
be easily extended to our case. Introducing symmetric and antisymmetric combinations of the bosonic phase fields, 
0jy s and 6js, (2^j24|] which represent the bosonic sound-like excitations of the system in the gapless phase, the 
low-energy Hamiltonian reads: 



H 



a=S,A 



47T 



dx 



-^(vQj, a ) 2 + K a (ve N . a ) 2 



/ dx {cob {Qif jA ) {A + B [cos (29 JiS ) + cos (2Qj iA )]} + C cos (29 j. s ) + D cos (29 JA )} , 
(nay J 

where a is a short distance cut-off, and A, B, C, and D are coupling dependent (but cut-off independent) constants. 

Let us consider the A-sector first. The A-sector can be gapless only if l/y/K A > 2 and, simultaneously, 2y/K A > 2. 
This is, of course, impossible. Thus, the A-sector flows to strong coupling and develops a gap. p3[ 

The Hamiltonian for the S-sector, renormalized by the A-sector, will be of the form: 



H S = 



VS_ 
47T 



dx 



1 



(ve,/, s ) + k s (ve^s) 



Vs 
(■Ka)'- 



J dx cos (29,7,5) 



The cosine term is relevant and opens up a gap when Kg < 1. Thus the system undergoes a KT transition when 
Ks — * 1. As discussed in the following sections, this is associated to a roughening transition. 

If the symmetric sector flows to the free Klein-Gordon Hamiltonian (i.e., the Luttinger model) in some range of 
parameters, the low energy spectrum of the two chains, for L — > 00, will have the form of a spinless Luttinger model 
describing symmetric excitations. Expressing 9j^^ s m terms of bosonic creation operators we can write: H] 



k 



(14) 



where Ns and J5 are the symmetric sector total number (N) and current (J) operators, vs is the renormalized sound 
velocity, v N = v s /K s , and vj = K s v s . @ 

In order to compute ujv, we note that the simplest charge excitation not involving the current part consists in 
adding 2 particles to the system. Thus we have: 

v N (L) = - [E(AN = 2) - E(0)] . 



To compute vj, notice that if a magnetic flux $ is concatenated with the ring, the current part of the energy spectrum 
is modified in this way: 



7T 

4L 



V J J S 



K ( T A ^ 

— vj Js + 4 — 
4L J \ b $ 



($0 is the elementary flux quantum). Therefore: 



vj(L) = 



L d 2 E(<$>) 
8^<9($/$ ) 2 



Finally, vs can be computed from 



v S (L) 



L 
27 



[E(k = 2tt/L) - E( 



where E(0) is the ground state energy and E(k — 2n/L) is the energy of the lowest excited state of momentum 
k = 2n/L. 

As a consequence, Ks can be equivalently computed from the finite-size extrapolation of vj(L)/vs(L), of 
vs(L)/vn(L), or of y/vj(L)/vN(L). If the finite-size data are compatible with a Luttinger Liquid picture, i.e., 
with a spectrum of the form (LL4|), then these three extrapolations should converge, as L — » 00, to a single value. 
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VI. ORDER PARAMETERS 



We now define the order parameters and correlation functions we have to consider in order to study the phase 
diagram of our model. We will basically deal with four correlation functions, whose behavior in the different phases 
is summarized in Tab. 1. 

• The height-height correlation, defined by Gh{r) — {{h r — ho) 2 }, diverges logarithmically as 

IK 

G h (r) = ln(r) + ... 

in the rough phase, with K > 1, while it remains limited in the flat, Neel, and DOF phases. At the 
roughening transition, K takes the universal value of 1 [p]]!^]. This gives a simple criterion for determining 
whether a phase is rough or not. In fact, the coefficient K coincides with the Luttinger exponent Ks for the 
symmetric sector, jl6| which can be extracted by finite-size scaling of exact diagonalization data (cfr. Sec. [v|). 

• The string correlation function, defined by (^] 

r-l 

G s (r) =-((S z exp(*7r]Ts;)^} . 



(We introduced the notation S? = — n^.) The phase factor contributes a plus (minus) sign if there are an 
even (odd) number of steps between site and site r. In the DOF and in the Neel phases, a step up (down) 
is preferentially followed by a step down (up). In these configurations, G s (r) gets a contribution equal to 1 
every time sites and r are occupied by a step. Thus, in the DOF phase and in the Neel phase G s (r) decays 
exponentially to the square of the mean density of steps. JT^ j 

• The staggered magnetization, defined by 

Af= lim I £(-!);<<?* <?*> . 

j 



A Neel phase will be signaled by a non-zero staggered magnetization M, while Af is zero in the rough, DOF and 
flat phases. 

• The flatness order parameter, defined by ||] 

r j=l 



T has a non zero value only in the ordered flat phase (in the DOF phase the exponential fluctuates between 1 
and —1 as r is increased). 





Flat 


Rough 


DOF 


Neel 


K s 


gapped 


> 1 


gapped 


gapped 


Gh (r) (for large r) 


< oo 


_» l n ( r ) 


< oo 


< oo 


G s (r) (for large r) 


-> 


-> 






j\f 


= 


= 


= 




T 


> 
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VII. THE OVERALL PHASE DIAGRAM 



Our model, even if the interactions are truncated to first neighbors, contains many parameters. Rather than trying 
to describe the phase diagram in an exhaustive form, we will focus our discussion on a few questions, that we consider 
quite relevant with respect to the surface physics interpretation of our model. This will take us to consider in detail 
some special planes cut through the phase diagram, and will also give us an idea of its global structure. 

Let us consider, once again, the Heisenberg spin-1 phase diagram in Fig. [|. From the surface physics point of view, 
it presents some unpleasant features: since the increasing temperature curve, for a given surface, is a line through the 
origin (the origin corresponds to the infinite T point; cfr. Eq. @)), every "surface" with repulsive interaction between 

steps of the same kind (V^ 11 > 0) has a preroughening transition at finite temperature, and no rough phase at finite 

T. On the other hand, if is attractive there is only roughening. 

In relation to these problems, we will discuss the following main questions: 

A Is the attractive Vf 1 term between opposite sign steps essential in order to stabilize a DOF phase? 

B Is there a choice of the parameters for which our model can describe a surface with a finite roughening temperature? 

C How does the presence or absence of a preroughening transition depend on the relative strength of the step-step 
interactions to the cost per unit length of a step. 

D What is the role of the opposite step on-site repulsion Vq~. 
In the following we address these points directly. 



A. The role of attraction between opposite steps: Spin-1 chain with V 1 — 0. 

In order to explore the different roles of the two interactions Vf 1 and vj, we have first studied the effect of a 

repulsive V}\ keeping Vi~ = 0. For the time being, we still work with the spin-1 condition, i.e., we impose an infinite 
on-site repulsion of opposite steps (Vq 1 = oo). 

For this choice of parameters, we do not find any point with V± > in which the finite-size data might indicate a 
vanishing gap. This is compatible with the results of den Nijs and Rommelse about the location of the KT-transition 
in the Heisenberg spin-1 phase diagram. Q The system undergoes a roughening transition only at infinite temperature. 

In Fig. |J we draw a qualitative phase diagram for values of the parameters \i and V" which are relevant for surface 
physics (i.e., positive energy cost for a step, and repulsive interaction between steps of the same kind). It is quite 
remarkable how the DOF phase survives the turning off of the attraction between steps of the same kind. As a matter 
of fact, taking V-y = leads to the disappearance of the reconstructed (Neel) phase from the physically interesting 
region of the phase diagram, and, therefore, to an even larger DOF phase. (We will further discuss the roles of 
and V]" in stabilizing the DOF phase later in this Section). 



B. Infinite or finite roughening temperature? The role of t ex . 

In order to discuss the point concerning the roughening temperature, we observe that, as T — > oo, the kinetic terms 
tend to be the only relevant pieces of the Hamiltonian, see Eq. (|J). Thus, we now consider the model in absence of 
potential terms (V), and for zero chemical potential fi. For the time being, we also take t ex — 0. The crucial role of 
the £ ea; -terni will be discussed afterwards. In this case, the Hamiltonian reduces to that of two coupled XY chains: 

H = + H.c) - 4 ;>>:., + H.c.) - ti ]T(S+A- la + H.c). (15) 

Notice that exchanging t" with is simply equivalent to renaming the sites (2i,a) to (2i,a) and vice- versa. This is 
illustrated pictorially in Fig. ||. Thus, a model with i^/t" = r is completely equivalent to one with i^/t" = 1/r. 
In Fig. ^(a) we plot the finite-size gaps as a function of the system size L for = 0, and different values of 
between and 1. Given the negligible curvature of the straight-line fits, the data seem to suggest that the gap 
extrapolates to in al l cases. In Fig. |?| we plot the finite-size value of the Luttinger exponent Ks, determined as 



explained in sec. VII B. The data for Ks confirm the scenario of a gapless (i.e., rough) system. Notice that K$ seems 
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to extrapolate to values larger than 1, indicating a rough phase that should survive to the turning on of a suitably 
small repulsive vj . 

We now address this point in more detail. Consider the = ill case (with t^ = 0), in which Kg seems to 



extrapolate to the largest value. Denoting by 
we can define the following four states: 



with s, s' = ±1, the four possible configurations at each site, 



IT> = 



II) = 



|o±) 



i 

71 



± 



It is now straightforward to verify that, at each site, the state |0_) is decoupled from the three remaining ones. The 
Hamiltonian H can then be considered a spin-1 Hamiltonian acting on the subspace spanned by the three states 
\]), | J,), and |0+). As it can be checked by explicitly calculating all matrix elements, H, when restricted within this 
subspace, coincides with the Heisenberg spin-1 Hamiltonian at J z = and J xy = 2t". As argued by den Nijs and 
Rommelse, Q the location of the KT transition in the Heisenberg spin-1 phase diagram is, very likely, exactly at 
J z = 0. Thus, two coupled XY chains with if- = t" and t$ = have a Ks actually equal to 1, and what we see in 
Fig. [t] is only due, very likely, to finite-size effects. 

The effect of turning on ij-, while keeping t^ = 0, leads to a completely different picture. In this case, for t^ = t" 
the system is gapped, as suggested by the finite-size data of Fig. ||(b). The physical reason for the different behavior 
of the t^ and ij} terms can be understood by considering the limiting cases of large values for these parameters. For 
to~ — > oo, the ground state tends to have |0+) at each site, with a large gap (of order ijj~) to other exited states. For 
tj; — > oo (at tjj- =0), on the other hand, the system reduces, by the previously described duality property, to two 
uncoupled XY chains, and must, therefore, be gapless. 

The previous considerations lead us to conjecture that for any choice of t^- (as long as t$ =0), two XY chains are 
gapless and have Ks = 1. On the contrary, turning on ig , at if" = 0, immediately opens up a gap. These conclusions 
have important consequences on the stability of the rough phase. Since Ks attains, at best, the marginal value of 1, 
turning on any positive vf immediately opens up a gap, and the rough phase is confined to infinite temperature. 

We will now demonstrate that, if we allow for the possibility of step-crossing events, t ex > 0, the gapless phase 
survives the turning-on of a positive vf , and every "surface" has a rough phase for high enough T. 

In order to show this, we add to the Heisenberg spin-1 Hamiltonian a t ex term, see Eqs. (f7|,|l2|), 

-t ex J2[(St) 2 (Si +1 f+ti.c.} , 



with t ex = ill = J xy . Fig. H shows the phase diagram for this case. Qualitatively, it is very similar to the Heisenberg 
spin-1 case (see Fig. H), except for small values of the potentials, where the t ex term changes the structure of the 

III 

phase diagram. In fact, for u = we observe a gapless phase extending for positive values of V 1 , up to V{ w 0.4: 
this is demonstrated in Fig. O were we plot the Luttinger exponent Ks along the line [i = 0. This finding is in accord 
with bosonization: the t ex term, unlike the t 1 - terms, increases Ks and leads to a stabilization of the rough phase. 
Indeed, for two XY-chains it is easy to show that, up to lowest order in t ex , 



K s = l + 

Another remarkable feature of the phase diagram in Fig. 
(t ex — 0), the temperature line for a given "surface" crosses t 

r\\ 



t, 



8ir 2 tU 



is that, at variance with the ordinary Heisenberg case 
ic DOF region only if the cost of a step, 5s, is sufficiently 

small as compared to the interaction between steps, . We have illustrated this by sketching in Fig. || temperature- 
lines for three different situations. For the case labeled A, the energy cost of a step is high with respect to the 
interaction energy between steps, and there is no preroughening. In the case labeled B, 5s /vf is smaller, and a DOF 
phase is present at intermediate temperatures. Finally, for case C, the interaction between steps is the most relevant 
energy, and the low temperature phase is 2 x 1 reconstructed. 



C. Presence or absence of preroughening: Role of interactions versus step line tension. 



We now want to discuss in some detail what happens if we leave the condition = —V^ , without going to 
the extreme case V-^~ = 0, discussed in Sec. VII A. We illustrate this by choosing Vf 1 = — V x /10, while keeping 
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t ex — i" =1 and = oo (infinite on-site repulsion of opposite steps). This choice of parameters describes a class of 
surfaces in which the attraction between steps of opposite kind is much smaller than the repulsion between steps of 
the same kind. 

In Fig. we plot the phase diagram for this choice of parameters. The system is Neel ordered for very large values 
of Vi, the value of the ratio Vi'/I^tH determines the location of the DOF-Neel phase boundary (see Fig. || and Fig. 
[|; recall that, for = 0, the Neel phase is absent for physical values of fi). The most relevant comment to the 
phase diagram in Fig. |l^ regards the conditions upon which the temperature trajectories of an actual surface model 
cross the preroughening line. It is clear, in fact, that depending on the ratio between the cost of a step (per unit 

length) 5s and, say, the interaction energy between steps of the same kind vf , a surface can have: i) only roughening 

(case A), or ii) first preroughening and then roughening (case B). Now, with = —Vf/10, the "critical ratio" 

(8s/Vi) CT it, below which preroughening is possible is of the order of 1/10, much smaller then in the Vf 1 = — Vj" case 

(where (#s/^i )crit ~ !)• Given the fact that 5s is typically the largest "diagonal" energy, this implies that a physical 
temperature trajectory will be, most likely, in the region where only roughening occurs. If, and how, long-range 
interactions might change this picture is an interesting and open problem. 



D. Role of opposite-step on-site repulsion: Finite Vq 



At last, we want to discuss briefly what happens to the DOF phase if we allow double occupation of a site, i.e., if we 
do not take the limit Vg 1 — ► oo. To demonstrate that the restriction to = oo is not crucial, we consider the case 

vf = — Vf 1 = 1, i" = ti = 1, with Vq~ finite. If Vp~ — > oo the system is an Heisenberg spin-1 chain at the isotropic 
point, corresponding to a DOF phase. || In Fig. [ll] we plot the finite-size values of the flatness order parameter T 
(open symbols) and of the DOF correlation function G s (L/2) (full symbols) for decreasing values of V^- . The data 
suggest that the system remains DOF all the way down to Vq ~ 0. We have verified that a similar scenario is found 
if we turn on the t ex term or a small Thus, our finite-size data suggest that the spin-1 condition (V^ — oo) is 
not essential in order to stabilize a DOF phase. 



VIII. STEP-STEP CORRELATIONS 



Correlation functions involving steps can be calculated numerically, for a given finite size, at any point in the phase 
diagram of our model. We will discuss here two correlation functions, i.e., step-step correlations and terrace width 
distributions. Let n s be the average density of steps of a single species (up or down). In general ns is always different 
from zero, even in the flat phase, since we do not discriminate between steps that traverse the entire sample and steps 
that form loops (i.e., finite terraces). Step-step correlations are defined as follows: 



7V T<T (r 



■(Step T (0)Step' T (r)) = — (n , T 7 



(16) 



with a =T, !■ If translational symmetry is not broken, we must have, at large distances, N^ a (r — > oo) — > 1. The 
distribution of terrace sizes (along the x-direction only!) is the probability of having two steps a distance r apart 
without any other step in between. There are two different kind of terraces we can look at: those delimited by two 
steps of the same type, and those between two different steps. Thus, we define: 



P T<T (r) 



1 



( n o, 



3 = 1 



(17) 



where, again, a =f, j- The string operator in square brackets enforces the absence of additional steps between and 
r. Fig. [l^ illustrates the behavior of and N^, at three different points in the phase diagram of the Heisenberg 
spin-1 chain: a rough case (J z = —0.5, /i = 0, triangles), a DOF case (J z = 1 and /i = 0, squares), and a flat one 
(J z = 1 and jU = —2, pentagons). The flat case results are very simple: both and N** converge exponentially fast 
(with a very short correlation length) to the large distance limit of 1. In the rough phase, instead, we have verified 
that the approach to 1 shows a power law tail. This is easy to prove. Rewrite first in terms of density and "spin" 
correlations: 



N^(r) = Jjg [(n n r ) ± (S^S z r 



(18) 
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where the + and — signs apply, respectively, to a =] and a =|, rii = n^j + ra^x, and Sf = Ui t -\ — n^- Within a 
bosonization approach, [^3|| the operators n, and involve (after particle-hole transformation for the |-bosons) only 
the antisymmetric and symmetric sectors, respectively. The antisymmetric sector is always gapped (see discussion in 
Sec. and Ref. HI), so that density-density correlations are exponential. In the rough phase, however, the symmetric 
sector is gapless, and S z — S z correlations have a uniform power-law tail of the form 

(SoSO = -T% + ••• , (19) 

which is precisely the term responsible for the logarithm in the heigth-heigth correlation function G(r) = ((h r — ho) 2 ). 
p5[ The DOF case results, finally, show a different behavior, with a sizeable oscillating component of the correlations. 
This behavior, however, reflects only a short-range effect, caused by the neighboring reconstructed (Neel) phase: the 
oscillating part has to decrease to zero at large r, since no breaking of translational symmetry occurs in the DOF 
phase. 

We finally discuss briefly the behavior of the distributions of terrace sizes (for simplicity, once again, in the Heisen- 
berg spin-1 case). While in principle it is important to know what is the probability for the surface to be flat over a 
distance r, this quantity has never been calculated so far. 

Let us consider, first, the behavior of P^l m the rough phase. Fig. |l3](a) is a plot of the logarithm of versus 
a scaled distance 2n$r, for several points taken inside the rough phase of the Heisenberg spin-1 phase diagram. We 
observe that the general behavior of P(r) is exponential in the size of the terrace r, 

P(r) pa er r/x 

and that a good collapse is obtained for all data if the distance r is scaled to the average separation between two steps, 
l/(2ns), i.e., A oc l/(2ns). The scattering of the data for the largest r's is due to finite-size effects. The behavior of 
P^(r) is found to be qualitatively similar. 

In the DOF phase we find that the terrace size distribution probability is again exponential with size, but now A 
does not scale with the density of steps, as it did instead in the rough phase. Fig. |l3|(b) illustrates the behavior of 
P{r) at a DOF point, corresponding to the isotropic Heisenberg point of the spin-1 chain. P(r) at a rough point 
is also reported for comparison. We observe that, as anticipated, the behavior of P'' 7 is, once again, exponential in 
r (i.e., the "DOF checkerboard" has no typical length!). Superimposed on the leading exponential, the DOF case 
results show a strong oscillating short range component which is again due to the neighboring reconstructed (Neel) 
phase. Two more features are worth noticing. First, compared to the rough case, P'^(r) is larger in the DOF case 
for r = 1, and then substantially smaller for larger values of r (and decreasing with a larger exponent). Second, in 
the DOF case P^ (r) is one order of magnitude smaller than P^(r), while the difference is much smaller in the rough 
case. These features are reasonable in view of the diluted antiferromagnetic ordering of steps, typical of the DOF 
phase. 

Experimentally, terrace sizze distributions could in the future be extracted, e.g., from STM data. |2q] 



IX. SUMMARY AND CONCLUSIONS 



In this paper we have presented and discussed a statistical mechanics model for studying the possible phase tran- 
sitions of an ideal, unreconstructed surface. The elementary objects upon which the model is based are the natural 
extended defects of an unreconstructed surface, i.e., steps and terraces. This starting point is, in our opinion, phys- 
ically more transparent then the usual microscopic RSOS-model description. Our model allows, in principle, the 
description of a real surface and, in perspective, one could test it with realistic step-step interactions. 

We have tackled our problem of interacting steps by mapping it, in a well known way, onto a one dimensional 
quantum problem of interacting hard-core bosons. Although this mapping is exact only in the strong anisotropy 
limit, it can provide very useful information about the phases and the nature of the transitions also in more general 
instances. Moreover, some realistic cases, like (110) surfaces of fee metals, are actually quite anisotropic. 

The quantum Hamiltonian, see Eq. (^), contains standard terms, like nearest-neighbor hopping (describing kinks 
on the steps), potential terms (describing interactions between steps), and chemical potential (cost per unit length of 
a step), as well as terms describing i) terrace creation/annihilation (through BCS-like number non conserving terms), 
ii) opposite step crossing events. The latter two terms are crucial, in many ways. 

Terrace terms are important to describe correctly the universality classes of the relevant transitions. This is known 
in the literature, p7| , ^6|jr^ ] but never explored in details in the present context. Moreover, in our case, the terrace 
terms also force us to work with hard-core bosons, as the standard Wigner- Jordan transformation to fcrmions does 
not lead to a simple local fermionic Hamiltonian. This point is sometimes overlooked in the literature. p8j 
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The term describing the crossing of opposite steps is important in order to stabilize a gapless (i.e., rough) phase for 
finite repulsive interactions between steps of the same kind. This, in turn, leads to a finite roughening temperature 
for the classical model. 

Finite-size exact diagonalizations and bosonization techniques have been used to unveil the richness of the phase 
diagram. In the limit of — > oo and for a particular choice of parameters (the potentials, for instance, are truncated 
to first neighbors and set to vf — — Vf 1 ), the model maps exactly onto the Heisenberg spin-1 chain Hamiltonian. The 
latter was also obtained, by den Nijs and Rommelse as the quantum mapping of RSOS models; || it presents a DOF 
phase for V} > 0, but does not describe, in that case, a surface with a finite temperature roughening. (On the other 
hand, if vf is attractive there is only roughening). 

Taking the Heisenberg chain as a starting point, we have then explored the phase diagram for other choices of 
parameters, obtaining results that we believe to be relevant with respect to the surface physics interpretation of our 
model. Summarizing, we have seen that: 

1. The Heisenberg spin-1 restriction v} = — is not crucial in stabilizing the DOF phase. In particular, we 



observe a DOF phase even for Vf~ = (see Sec. VII A ). Moreover, a DOF phase is present not only for V" = oo 



(spin-1 case) but also when is finite, as long as positive. 

If we add to the Heisenberg spin-1 Hamiltonian a t ex term, we observe a gapless phase extending for positive 
values of v} . Ever y surfa ce has a rough phase for high enough T. This is true also for other choices of the 



potentials (see Sec. VII C). Moreover, if we do not include in the Hamiltonian the t ex term, the rough phase 



does not survive when one turns on a V-l > (see sec. VII B). Thus, the opposite step crossing term is crucial 



in order to obtain a model describing, at least at a coarse grained level, a physical surface. 

3. The relative values of the interactions and of the cost per unit length of a step decide whether a surface has 
a stable DOF phase for a certain range of T. The temperature trajectory crosses the DOF region only if the 
cost of a step, 5s, is sufficiently small as compared to v} (see Fig. ||). Given the fact that 8s is typically the 
largest "diagonal" energy, this implies that a physical temperature trajectory will often be in a region where 
only roughening occurs. 

In conclusion, we have found that: (i) a model based on steps can describe prcroughcning (PR), as well as rough- 
ening; (ii) the steps must be treated as hard-core bosons rather than fermions; (hi) the qualitative role of step-step 
interactions in driving PR, known already from RSOS models, is recovered in this picture; (iv) correlation functions 
involving steps can be calculated in a quite straightforward way. The main one, never studied so far, which we have 
considered, is the terrace size distribution. Here we find simply an exponentially decreasing probability for increasing 
size. This result should be amenable to experimental testing, for example by STM; (v) In view of the additional 
simplicity of step models, it should be feasible, in the future, to study the role of long-ranged interactions, a problem 
without hope of solution within RSOS models. 
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FIGURE CAPTIONS 

FIG. 1. Scheme of a surface with up (j) and down (J.) steps. The heights of the terraces are explicitly indicated. The other 
symbols refer to the Boltzmann weights considered: 8k (cost of a kink), St (cost of a terrace creation), S ex (cost of a step 
crossing), V" and V 1 ' (interactions between parallel and opposite steps); the black dots indicate where terraces are created or 
destroyed. 

FIG. 2. Schematic representation of a kink, the beginning of a size-1 terrace, a size-0 terrace, and a step crossing between 
strip j and strip j + 1 and the relative energetic costs. 

FIG. 3. Phase diagram for the spin-1 Heisenberg chain. || Here and in the following we consider only negative values of 
the chemical potential fi, that are relevant for the surface physics problem. Lines (a) and (b) are discussed in the text. 



FIG. 4. Qualitative phase diagram for a spin-1 chain with t\ = t" = 1, t ex = 0, V\ = 0. 



FIG. 5. The duality mapping for two XY chains, see Eq. ([15j). Dotted, solid and dashed lines denote, respectively, t x , t" , 
and to couplings. The duality survives also in presence of Vq~ . 
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FIG. 6. (a) Finite-size charge excitation gaps for the two XY chains in Eq. ( |l5| ) at tg = 0, t" = 1, and for various values 
of t x . Dashed lines are obtained from straight lines constructed so as to pass through the L = 6 and L — 8 points. The 
extrapolation to zero is remarkably good, (b) Finite-size charge excitation gaps for i x = 0, i x = t" = 1 (dashed line), and 
t x = 0, to = t" = 1 (solid line). The solid and dashed straight lines are constructed as in (a). Notice the remarkable smallness 
of size corrections. 



FIG. 7. Luttinger exponent Ks for the same parameters of Fig. |^(a). As argued in the text (see sec. VII B), Ks should 



converge to 1 for L — » oo and the apparent extrapolation to values larger than 1 is very likely due to finite size effects. 

FIG. 8. Qualitative phase diagram for a Heisenberg spin-1 chain with exchange term, i.e., t x = t" = t ex = 1, = — V-^ . 
Lines (A), (B), and (C) are discussed in the text. 

FIG. 9. Finite-size Luttinger exponent Ks along the line fj, = 0, for the Heisenberg spin-1 chain plus exchange term. Ks 
extrapolates to values larger than 1 for Vj" < 0.4. Dashed lines are only guides to the eye. 

FIG. 10. Qualitative phase diagram for a spin-1 chain with i x = t" = t ex = 1, Vj^ = —Vj'/lO. Line (A) describes a situation 
where the cost of a step 5s is larger than the interactions, and only roughening is found. Line (B) describes a situation where 
8s is smaller than \V\ |, and roughening is preceded by PR. 

FIG. 11. Finite-size values of the flatness order parameter T (open symbols) and of the DOF correlation function G s (L/2) 
(full symbols) at the Heisenberg isotropic point, for decreasing values of V . The system appears to be DOF for all positive 
values of V X . The dashed lines are only guides to the eye. 

FIG. 12. iV n (ab ove) and (below), at three different points in the phase diagram of the Heisenberg spin-1 chain: a 
rough case ( J z = V ± = —V{ = —0.5, D — —fj, — 0, triangles), a DOF case ( J z = V ± = — V{ = 1, and D = —fj, — 0, squares), 
and a flat one (J z — V ± = —V} = 1, and D = —fi = 2, pentagons). Lines are only guides to the eye. 

FIG. 13. (a) In(P^) versus the scaled distance 2nsr at various points in the rough phase of the Heisenberg spin-1 chain; full 
symbols: J z = V ± = -V} = -0.25, D = —jj, = 0,0.1,0.2,0.3; empty symbols: J z = V 1 " = -V^ 1 = -0.5, D = —fj, = 0,0.2,0.4; 
stars: J z = V ± = —V} = —0.75, D = —fj, = 0,0.1,0.2,0.3. (b) and , in logarithmic scale, at two different points in 
the phase diagram of the Heisenberg spin-1 chain: a DOF case (J z — V 1 - = —V± = 1, and D — —fj, — 0, squares), and a rough 
case (J z = V ± — — V]" = —0.5, D = — fj, = 0, triangles). Full and empty symbols correspond to P Ti and P^ , respectively. 
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